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Abstract 

In this work we study Podolsky electromagnetism in thermodynamic equilibrium. We show 
that a Podolsky mass-dependent modification to the Stefan-Boltzmann law is induced and we use 



experimental data to limit the possible values for this free parameter. 
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I. INTRODUCTION 



Electromagnetism is a U(l) gauge theory lj. The electromagnetic field emerges from 
the application of the gauge principle to the local Abelian group. This principle tells us 
that the electromagnetic field must be a Lorentz vector with some well-defined internal 
properties. What the principle does not tell us is which is the highest order of the field 
derivatives appearing in the Lagrangian of the theory. Usually, we just invoke Occam's razor 
and use the simplest one: a Lagrangian with only first-order derivatives. It is well known 
that there exist Lagrangian densities of second-order derivatives which are equivalent to 
Lagrangian densities with derivatives of only first order. For instance, consider the following 
two Lagrangian densities for the real scalar field: C\ = ^d^&d^ + V^) and £2 = — + 
V(<f>). They both lead to the same equation of motion [2]. Since they lead to the very 
same physical consequences, we are free to choose either one or the other. This situation 
is unlike what occurs with the electromagnetic field. It has been proved that there is only 
one extension (up to a total-divergence term) for a second-order derivative Lagrangian for 
electromagnetism that is both Lorentz and gauge-invariant This extension is known as 
Podolsky electromagnetism, after Boris Podolsky has first proposed it in 1942 [4]. As we shall 
review in the next section, Podolsky Lagrangian does not lead to the equations of motion 
expected from the Maxwell theory. Therefore, they are non- equivalent descriptions of the 
Abelian gauge field. 

Podolsky electromagnetism depends on a free parameter. The way used to fix the Podol- 
sky parameter is the same one used to set bounds on the values of all free parameters of 
the Standard Model of Elementary Particles: it can only be fixed from experiments jsj]. In- 
deed, this is one of the aims of this paper. Besides, despite its long-dated success, Maxwell 
electromagnetism still has unsolved problems. In the classical level, for example, the elec- 
trostatic potential diverges over punctual electric charges. On the other hand, Podolsky's 



electrostatic potential is finite everywhere [4J. In addition, using Podolsky theory, Frenkel 
was able to solve the famous "4/3 problem of classical electrodynamics" jf|, while in the 
context of Maxwell theory this problem remains open. Furthermore, Podolsky theory of 
electromagnetic interaction presents a richer theoretical structure than a theory with only 
first-order derivatives. Finally, with Podolsky electromagnetism there is the possibility of 
new Physics. We see, then, that Podolsky theory is interesting from both theoretical and 
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experimental points of view. 

In 1901, Planck's pioneer and now famous work on black body radiation led to the 
foundations of Quantum Physics. In order to fit the available data, Planck postulated 
that the energy exchanging between cavity oscillators and the Maxwell electromagnetic field 
in thermal equilibrium is quantized. Under these assumptions he derived the frequency 
distribution of the black body radiation. In his work, Planck was also able to deduce the 
Stefan-Boltzmann law, which states that the power for unit area of the black body radiation 
grows with the fourth power of the temperature. We can cite as an example of Planckian 
black body distribution the Cosmic Microwave Background Radiation (CMBR) 7]. In fact, 
with precise data, CMBR was considered the most accurate black body radiation measured 
up to date 8|. Nowadays, Planck's law can be derived using modern methods of Quantum 



Field Theory at Finite Temperature 



ll| . In this context, black body radiation is seen 



as a gas of Maxwell photons in thermodynamic equilibrium. In the light of this thought, in 
the present work we investigate a gas of Podolsky photons at thermal equilibrium and we 
seek the finite-temperature properties of the Podolsky theory. Through use of imaginary- 
time technique from Finite- Temperature Field Theory we construct the partition function 
of the theory after determining its constraint structure. From the partition function we are 
able to evaluate all thermodynamical quantities, including the energy density distribution. 
We expect a modification on the Stefan-Boltzmann law due to the presence of the term that 
contains higher-order derivatives in the Podolsky Lagrangian density. Finaly, we compare 
our results with experimental data for the Stefan-Boltzmann law at CMBR temperature 
and we set a thermodynamical limit on the Podolsky parameter. 

This paper is organized as follows. In section [III we review the canonical structure of 
Podolsky electromagnetism. In section II III we work on the path integral formulation of the 
transition amplitude which is a step to the attainment of the Partition Function of the 
problem. In section IIVI we use imaginary-time formalism and deal with the problem of 
evaluating the Partition Function. In that section we also show the Podolsky correction to 
the Stefan-Boltzmann law and we use experimental data to limit the free parameter of the 
theory. Our final remarks are presented in section [V] 
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II. GENERAL ASPECTS OF PODOLSKY THEORY'S CANONICAL STRUC- 
TURE 



The content of this section is by no means new. We intend only to review some general 
aspects of the canonical structure of Podolsky electromagnet ism. 
The Lagrangian density for Podolsky theory is 1 4j 

C = ~F lw F ,w + J-d u F^d,F^ 



(1) 



dyA^. The free 



' 2m 2 M 

where, like in Maxwell's case, the field-strength tensor is F^ = d^A 
parameter m has dimension of energy and in the limit \m\ — > oo we recover the Maxwell 
Lagrangian density. The Euler-Lagrange equations for this theory are 



(□ + m 2 ) dpF*" = 0. 



(2) 



As we have anticipated in the introductory section, these equations of motion differ from 
those of Maxwell theory. Therefore, the Physical contents of the Podolsky theory do not 
coincide with those of Maxwell theory. 



The symmetric Energy-Momentum density tensor reads 
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1 



(2d x F^d x F» - 2d x F^d i F x v + d x F x ^d^ u ) 



(3) 



The energy density S is the component 00 of this tensor. It is possible to write £ in terms 
of the electric and magnetic fields: 



£ 



1 



E + W 



(V-E) z + E-VxB +4E-QE + 4E- V (V • E) 



(4) 



2 l m z 

This expression does not appear to be positive-definite in the general case. However, if 
we restrict it to the electrostatic case, we have for the energy [4j: 



E, 



electrost. 



d 3 x£ 



electrost. 



(Fx 



E 2 + \ (V ■ E) 2 

m 2 



(5) 



1 Throughout this work we set natural unit system, the metric signature (- 
as we work in the Minskowski spacetime. 



-), and □ = d^dfj, as long 
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Once we impose the condition E e i ectTOSt . > 0, we have the implication that the parameter 
m must be real and, without loss of generality, we assume it to be positive. 

Now, we can impose the generalized Lorenz condition (□ + m?) d^A^ = on the Podolsky 
field and the equations of motion ([2]) become simplified ljj: 



(□ + m 2 ) OA* = 0. 



(6) 



One possible solution to this equation can be written as 



A» = A» M + A P , (7) 

where A^ M satisfy the equations of motion of Maxwell electromagnetism OA^ t = and A p 
satisfies the Proca equations (□ + m 2 ) A p = 0. Here we notice a notable difference between 
Maxwell and Podolsky theories. Although Maxwell field has only four components (one for 
each spacetime direction), Podolsky's has eight. For each direction the Podolsky vector field 
has one massless and one massive sectors, as it is seen in (J7J). This interpretation will be 
useful in the next sections. 

We can write the canonical Hamiltonian of Podolsky theory He as 



H c = I d 3 x(p,A» + n,A»-£), (8) 
with the momenta defined as 



2V 



■n,. 



dC 



d [A» 
dC 

d(A» 



do 


dC 


- d k 




d (lA 





dC 



+ 



dC 



d(d d k A») d ( dk A» 



(9) 
(10) 



which are canonically conjugated respectively to the fields A^ 1 and y? M = A 11 , which is con- 
sidered as an independent variable. 

Expliciting these quantities [He = J d 3 xHc) '■ 
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P, = - Fop + (rffrdjF* - d d u F^) ; (11) 



1 



- = - r)\d k F\) ; (12) 



A< 9 



^ c = - p ■ ^ - im 2 ^ 2 + Tx k djF^ k - ~ + VA ) 2 + 

' (V 2 A° + V- V) 2 + -^' fc ^ fc - (13) 



2m 2 v ' I 

Now, we can follow the steps of reference 



13( for the constraint analysis a la Dirac 



and find that there are three (unlike only two found in Maxwell theory) first class constraints: 

^sttowO; (14) 
<j) 2 =p -V (15) 
3 =V-p^O, (16) 

with the usual notation the symbol "~" means weak equality. Following Dirac's procedure, 
we choose three gauge conditions: 

fii^^wO; (17) 
tt 2 = (D + m 2 ) V- A^O; (18) 
Q 3 =A°^0. (19) 



It was shown in [l3j that this set constitutes an appropriated non-covariant gauge con- 
dition which fixes the first-class constraints. Solution (J7]) has made it clear that Podolsky 
electromagnetic field has eight apparently independent components. However, gauge condi- 
tions (I17H19I) reduces this number to five. Five "degrees of freedom" (d. o. f.) are compatible 
with the interpretation that Podolsky field is composed of a Maxwell field (with two d. o. 
f.) plus a Proca field (which has three d. o. f.). 

III. PATH INTEGRAL FORMALISM 



Once constraint analysis is done, we are able to write down the generating functional 
with null sources (or transition amplitude): 
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Z = DprDTTsDA'Dp* det {Q a , b } 



n=l 



exp ( i I d 4 xC c ) (20) 



where 



det {tt a , <p b } = det [(□ + m 2 ) V 2 ] ; 

Cc=P»dtA» + -He. 



(21) 
(22) 



After some steps not quite different from those of Maxwell's case, we can rewrite the 
transition amplitude as 

Z = I Yl DA ° det [( n + ™ 2 ) V 2 ] 5 [(□ + m 2 ) V • A] exp fi f d 4 xA (23) 

with C given by equation p]). 

Using a straight-forward generalization of the Faddeev-Popov ansatz, we can pass from 
a non-covariant gauge fixing to a covariant one: 



.cr=0 



det 



(□ + m 2 ) □ 



(□ + m 2 ) 9 5 A" - f 



exp « 



d A xC 



(24) 



where p 7^ is an arbitrary real number and / = f(x) is an arbitrary real function. Physical 
quantities are independent of the function f(x). So, we "sum" Z over all functions /, 
considering the weight factor exp \—\ J d 4 xf 2 ): 



DfZexp(-- I d 4 xf 



I 



DcDc 



Y[da c 



cr=0 



exp I i I d 4 xC 



eff 7 



(25) 



where c and c are ghost fields. The effective Lagrangian density £ e ff is defined as 



' d^d^Fi - -L [(□ + m 2 ) d^] 2 - -c (□ + m 2 ) Uc. (26) 



2m 2 



2p 2 



P 



So far we have studied Podolsky theory at zero temperature. In the next section we will 
analyze the theory in thermodynamic equilibrium. 



IV. FINITE TEMPERATURE 



It is possible to obtain the Partition Function from the transition amplitude (1251) . Once 
the Partition Function is carried out, all thermodynamical properties of the system becomes 
available. In order to obtain the Partition Function from the transition amplitude we just 
have to perform an Euclideanization of the time components of the vector fields, a compact- 
ification of the Wick-rotated time coordinate and impose periodic boundary conditions (P) 
in this coordinate for the electromagnetic and the ghost fields lp|. 2 Doing so, we find the 
Partition Function for the free Podolsky field: 



Z{0) = I DvDc 



Y[ DA ° 



(7=0 



exp ( — J dxC,E 



(27) 



where /3 = 1/T, T is the temperature, and we use the notation 



dx = dr d 3 x; 

(3 Jo J 

A = — <9 7 <9 7 ; 

1 1 

C E = - -F CX F (X - — -<9 c F fA c\F 5A + 



2p 



1 [ ( A + m 2 ) d x A x ] 2 - -c (A + m 2 ) Ac. 



P 



(28) 
(29) 

(30) 



Ce is called the effective Euclidean Lagrangian density. Since there is no coupling among 
the electromagnetic field and the ghost fields, the Partition Function takes the form 



x 



DcDc exp < — I dx 

I J/3 

f[DA 

(7=0 



(A + m 2 ) At 



exp ( — - J dxA a M m A 1 



det 



P 



(A + m 2 ) A 



[Det {M ai )\ 



(31) 



2 The Euclideanization follows as a; — * —it (and therefore da —> id/dr = id T ) and A — * —iA^ = — iAq. 
All summations are performed with Euclidean metric from now on 15 1. 
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where "det" stands for the determinant in the Hilbert space alone (as usual) and "Det" 
stands for the determinant in both Euclidean spacetime and the Hilbert space. The operator 
M„~, is defined as 



M ai = I 1 



A 



m 



A 



1 + 



/ 2\ 2 



V P 



d a d^. 



(32) 



After evaluating the determinant in the Euclidean spacetime we find 



Z(p) = [det (A)]" 1 [det (A + m 2 )] 2 . (33) 

We note that the Partition Function is a product of determinants of the form 

-Hi 

[det (A + m 2 )] 2 , with j = 1, and 2. Each of these terms describes a gas of free par- 
ticles with mass rrij and rij d. o. f. We identify the first of these determinants as a partition 
function for massless particles with two d. o. f., i. e., Maxwell photons. On the other 
hand, the second determinant is the partition function for particles of mass m and three 
d. o. f. Those are Proca particles. Since Z ((3) involves no other terms, it describes a gas 
formed of non-interacting gases of free Maxwell photons and free Proca bosons. The "non- 
interaction" property between the two distinct gases is a direct consequence of Podolsky 
theory's linearity. 

In order to evaluate the determinants we note that the equation 



det (A + m 2 ) = J] (3 2 [u 2 n + u? (p, m)] , (34) 



n,p 



with the relativistic energy- momentum relation uj (p, m) = a/p 2 + m 2 and the bosonic Mai- 
subara frequencies u n = 2mr//3, remains valid for both massive and massless cases 11]. 
Using this identity, the logarithm of the Partition Function reads 



In [Z (P)] = - ]T In [/3 2 K + p 2 )] - 2 J> {/3 2 [ul + ^ (p, m)\ } . (35) 

n.p ra,p 

After evaluating the sum in n, passing to continuous in momentum space, and discharging 
irrelevant /3-independent terms and vacuum contributions, we have: 
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In [Z 09, V)] =-2V [Ms In (l - e"*) - 3V / In [l - e*-^] . (36) 

The first term in the r. h. s. of ( 1361) is associated with a gas of free Maxwell photons 
and it gives the Planck's law. Since it can be found in many text books we shall skip its 
computation. The second term, on the other hand, corresponds to the massive sector of the 
theory. This term depends on the Podolsky parameter and shall give a correction to the 
Stefan-Boltzmann law. Calling the m-dependent term as In (Z') we find, after changing the 
integration variable to x = uj/m: 3 



ln(Z' 



m A (3V 



oo 

E 

k=l 



dx(x 2 -l) 2 X * e -W™. 



(37) 



In order to write this expression in a compact form we use the following representation 
of the Modified Bessel Function of the Second Kind, valid for n > —1/2 161 ]: 



7T (Z\ n l°° „ , o .xn-4 



K n (z) = r £ + g (- j jf e—^-lJ^'dx, (38) 
where r(?/) is the Gamma Function. We can now write 

, jvv^w (39) 

As far as we know, there is no known analytical, closed form to the summation appearing 
in this equation. Hence, we restrict ourselves to evaluate In (Z') only approximately. We 
recall that in the limit m —>■ oo the results found using Podolsky electromagnetism must go to 
the usual results of Maxwell theory. For this reason, we expect In Z' to be a correction to the 
Planck law. In thermodynamic equilibrium, temperature is a natural scale of energy. In this 
sense, we will evaluate the correction to Stefan-Boltzmann law in the regime f3m — m/T ^> 1. 
So, we can write 



K'^-Jwk*'*" (40) 



The sum over k arises from the expansion of (l — e ,3w ) 1 . 
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and keep only the first term of the sum in (139]) . Within this approximation, we can solve 
equation fl36l) : 

TOV;™)] = ^ + 3^)V™ (41) 

The first term in the r. h. s. of this equation is the usual Planck result. The second 
term is a correction due to the Podolsky mass. This equation enables us to evaluate any 
thermodynamical quantity. The energy density u (T; m), for instance, is found to be 



ld\n[Z((3,V;m)] 

u (T; m) = — 



= a (T, m) T 4 , (42) 

v 



V dp 

with a (T, m) = o~q + 5a (m/T), where o"o = 7r 2 /15 is the Stefan-Boltzmann constant and 

^S=^ra f e- (43) 



TJ V8t^\T. 

is the correction due to the Podolsky parameter. As we can see, in the limit m/T — > oo, 
a (T, m) goes to o"o and we recover the Stefan-Boltzmann law in (]42l as expected. 

Equation ( f43l) can be used to set a limit on the possible values of Podol- 
sky parameter. The experimental value for Stefan-Boltzmann constant is o"o = 
(5,670277968 x 10" 8 ± 4 x 10" 13 ) W/(m 2 K 4 ) [17]. Since so far we have not detected any 
sensible deviation from the Stefan-Boltzmann law, the correction 5a (m/T) must be at most 
equal to the experimental error in the Stefan-Boltzmann constant. Our results showed that 
such a correction depends on the temperature of the black body radiation. Recalling that 
CMBR shows a very accurate black body spectrum, we set T = 2, 725K (the temperature 
of CMBR) and we find that all values for the Podolsky mass such that m > 4, meV are 
compatible with the experimental data. 



V. FINAL REMARKS 



In this work we studied the Podolsky theory for electromagnetism at finite temperature. 
We have reviewed the canonical structure of the theory and have noted that even though 
both Maxwell and Podolsky Lagrangians are possible choices for the Abelian gauge group, 
they give rise different physical results. Using imaginary-time formalism we showed, through 
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Partition Function evaluation, that a gas of free Podolsky photons in thermodynamic equi- 
librium is mathematically equivalent to a gas formed of non-interacting gases of free Maxwell 
photons and free Proca bosons in thermal equilibrium as well. We argue that this result is 
only valid for temperatures well bellow the one that corresponds to the electron rest energy 
If the temperature is raised enough, we should take into account fermion pair creation 18]. 
Of course, this can only be accomplished in Podolsky QED. We also showed that Podolsky 
electromagnetism induces a modification in Stefan-Boltzmann Law. Accordingly to Podol- 
sky theory at finite temperature, the existing corrections to the Stefan-Boltzmann law vanish 
as the Podolsky mass goes to infinity recovering, in this approximation, the original law as 
expected. Using experimental data for Stefan-Boltzmann constant and the temperature of 
the CMBR we set a thermo dynamical limit to the Podolsky parameter. The reason why we 
have chosen CMBR temperature is twofold. First of all, CMBR has been called the most 
accurate black body radiation. Second, it is a temperature well below that associated with 
the electron mass. Therefore, our results are meaningful in that regime. Our analysis have 
shown that the Podolsky mass cannot be smaller than approximately 4, meV, otherwise it 
would already have been detected in black body radiation experiments. We end this section 
stating that Podolsky theory remains as a possible choice for the electromagnetic field and 
only further work in both theoretical and experimental research fields will be able either to 
confirm it or to rule it out. 
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